Quantum anharmonic oscillators: a new approach 
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Abstract. The determination of the eigenenergies of a quantum anharmonic 
oscillator consists merely in finding the zeros of a function of the energy, namely the 
Wronskian of two solutions of the Schròdinger equation which are regular respectively 
at the origin and at infinity. We show in this paper how to evaluate that Wronskian 
exactly, except for numerical rounding errors. The procedure is illustrated by 
application to the gx 2 + x 2N (N a positive integer) oscillator. 
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Quantum anharmonic oscillators have been frequently used in different branches of 
Physics to simulate a great variety of situations and to explain multitude of phenomena. 
Apart from this, since the publication of the seminal papers by Bender and Wu PQ and 
by Simón and Dicke j2] showing the failure of the Rayleigh-Schròdinger perturbation 
method, they have served to test plenty of approximate methods of solution of the 
Schròdinger equation. Papers dealing with the most recently proposed methods jÜ] 
contain references to older ones, that we omit for brevity. It seems, however, to have 
been passed unnoticed that an exact procedure exists to obtain a quantization condition 
that gives the eigenenergies as zeros of an easily calculable function. The idea is the same 
one exploited in the solution of the Schròdinger equation for the harmònic oscillator or 
the Coulomb potential, although the procedure is slightly more tricky than in those two 
simple cases. 

The starting point is the Schròdinger equation written as a second order differential 
equation free of first-derivative terms. Such equation is satisfied by the wave function, 
in the case of an even one-dimensional anharmonic oscillator, or by the reduced wave 
function, if one is considering an isotropic D-dimensional oscillator. A solution u?eg of 
the differential equation physically acceptable at the origin can be immediately obtained 
as a power series of the variable. Two other solutions u^ 1 ' and vS 2 \ characterized by their 
behavior at large distances, can also be considered. To be specific, let u^ 1 ' represent the 
solution going exponentially to zero as the variable increases, whereas corresponds 
to an exponentially diverging one. Since vS^ and are independent, «reg can always 
be written as a linear combination of them with coefficients, called connection factors, 
which depend on the energy and on the parameters of the potential. For a genèric value 
of the energy, both connection factors are different from zero and «reg is not a physical 
solution because of its behavior at infmity. The eigenenergies are then determined by 
requiring the cancellation of the connection factor multiplying u^. 

The connection problem for a differential equation with two singular points (let 
us say, one at the origin and the other at infmity) was discussed by Naundorf 0j. He 
considered the case of one of the singular points (that at infmity, for instance) being 
irregular of integer rank R > and the other one being either irregular of integer 
rank r > or regular (r = 0). Here we are interested in this last case. Naundorf 
gave a procedure consisting in obtaining 2R independent formal power series, with an 
integer index running from — oo to +oo, having well defined asymptotic behaviors, and 
whose coefficients can serve as a basis in the 2i?-dimensional space of solutions of the 
recurrence obeyed by the coefficients of u Teg . To obtain such basis, Naundorf replaces, 
in the known expressions of and [5], the exponential term determining their 
respective asymptotic behaviors by R independent formal expansions of the type of the 
Heaviside's exponential series. Multiplication of those formal expansions by the Taylor 
series of the rest of the exponential terms and the descending power series in vP*' and 
produces 2R formal expansions whose coefficients obey the above mentioned recurrence, 
i. e., the required basis. Comparison of 2R consecutive coefficients of the power series 
expression of u reg with the analogous coefficients of the elements of the basis leads to 
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a system of 2R linear equations whose solution allows one to obtain the connection 
factors. That procedure has been applied to the solution of several physical problems, 
like the hydrogen atom with fine structure [Oj, the quarkonium [Jj, the spherical Stark 
effect in the hydrogen atom jH] or the quartic and sextic anharmonic oscillators jÜ]. 

The method here suggested is related to the Naundorf's one insofar as it also 
rests on the vanishing of one of the connection factors and makes use of Heaviside's 
exponential series to obtain formal expansions, but differs from the Naundorf's method 
in the procedure of computation: instead of following the steps detailed in the preceding 
paragraph, we benefit from the fact that the connection factor multiplying is given 
by the quotient of Wronskians W[u re g, u^/Wlu^, u^] and, since the denominator 
does not vanish, the quantization results from the fulfillment of the condition 

%c g ,« (1) ] = 0. (1) 

To implement this condition, we need suitable expressions of u Teg and The series 
expansion mentioned above is adequate to represent the first of these solutions. For 
the second one a closed expression does not exist, in general, but a formal (asymptotic) 
expansion can be easily obtained by substitution in the differential equation. Then it is 
trivial to write a formal expression of the Wronskian and to require its cancellation. 

To illustrate the method, let us apply it to the determination of the eigenenergies 
of the one-dimensional anharmonic oscillator represented by the potential 

V(x) = g x 2 + x 2N , iV a positive integer. (2) 

This problem has been tackled by several authors ^01 EU H2j by using different 
approximations. We discard the trivial case N — 1. The cases N = 2 (usually referred 
to as quartic oscillator) and N = 3 (sextic) can be easily solved following the steps 
we are going to detail, but the resulting equations do not fit in the general form given 
below. Therefore, we assume N > 4. The Schròdinger equation (in adequate units for 
the variable x and the energy E) 

admits solutions, regular at the origin, of the form 



1t rC g yX J 



£a n x™ + ", a ^0, (4) 

n=0 

with v = (even states) or 1 (odd states). Alternatively, two independent solutions, 
and u^- 2 \ with asymptotic expansions (for x — > +oo) 

" WI 1 -" , t^-<- "'■ 1<T / J = 1,2, (5) 



u^\x) ~ exp 



x N+1 



m=0 



N+l 

can also be considered. Substitution of this formal expansion in (j3J) gives for the 
exponents 

= -1, =H= -N/2, 

a (2) = +i ; ^(2) =fi = _ A r/ 2 , [ü) 
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and for the coeficients 

2aG>m/# = (m-N/2)(m-N/2-l) h%_ N _ x + Eh^ N+1 - gh { £_ N+3 .(7) 

Instead of computing directly the left hand side of (JTJ), let us introduce two auxiliary 
functions 

v mg (x) = exp (x N+1 /(N + 1)) u Icg (x), (8) 
v {1 \x) = exp (x N+1 /(N + 1)) (9) 
which obey the differential equation 

*1 - 2 x » ^ + ( E - g x 2 - Nx"- 1 ) v = 0, (10) 
dx 2 dx v ' 

and whose Wronskian satisfies 

WK^' 11 ] = exp (2x N+1 /(N+ 1)) WK s ,m (1) ]. (11) 
Now, by using the series expansion 

oo 

v reg {x) = Y J b n x n+u , b ^0, (12) 

n=0 

with coeficients given by the recurrence 

(n + v)(n + v-l)b n = -E 6 n _ 2 + g ò n _ 4 + 2(n- N/2-l + u) 6 n _jv_i , (13) 
and the asymptotic expansion 

oo 

~ E ^" m+M > ( 14 ) 
m=0 

one obtains for the left hand side of (jlljl a formal expansion 

oo 

W[« rcg yV E T*** -1 "^", (15) 

£:=— oo 

with coeficients 

oo 

7fc = E (~ 2m -k-v + fi) b m+k h m . (16) 

m=0 

A similar expansion can be obtained for the right hand side of (jllj) by recalling the 
Heaviside's exponential series 

oo t nJr ^ 

introduced by Heaviside in the second volume of his Electromagnètic theory (London, 
1899) and probed by Barnes P3j to be an asymptotic expansion for arbitrary 5 and 
| arg(í)| < 7r. Let us construct N + 1 expansions 

" (2x N+1 /(N + l)) n+s ^ . . 



ex V (2x N+1 /(N+l)) ~S L = Yl 



T(n + 1 + 5 L ) K ' 

of the type (fTTjl with appropriate choices for 5, 

5 L = {u + fi + L)/(N + l), L = 0, 1, . . . , N. (19) 
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Table 1. The four lowest eigenenergies of the oscillator for N = 4 and several 
vàlues of g. The energies Eq and E 2 correspond to even states (u = in (@J); i?i and 
-E3 to odd ones (y = 1). 



g E Ei E 2 E 3 



-20 


- 15.62781790 


- 15.60342843 


- 1.99759674 


0.04913769 


-10 


-3.89894214 


- 3.32541335 


3.26415045 


8.82212629 


-1 


0.93527862 


4.11346827 


9.49008984 


16.49163253 


-0.1 


1.19798114 


4.69299658 


10.16968229 


17.25807961 





1.22582011 


4.75587441 


10.24494698 


17.34308797 


0.1 


1.25340643 


4.81845727 


10.32015025 


17.42806187 


1 


1.49101990 


5.36877806 


10.99373734 


18.19110002 


10 


3.21296474 


9.86889192 


17.20002166 


25.52311499 


20 


4.48741520 


13.54543209 


22.89430780 


32.78247104 



It is evident that, for any set of constants /3l (L — 0, 1, . . . , N) restricted by the condition 

N 

= £/?l, (20) 

L=0 

one has 

N 

exp (2x N+1 /(N + 1)) W[u ies , u«] ~ J2 fa S l . (21) 

L=0 

If, according to Eq. (fTTJl . this formal expansion has to coincide with that in (|I3j). the 

constants (3i must be 

A = SInIvT* lk "' "l = ^+D + 1 + L, (22) 

where the integer n can be chosen at will. Substitution of those vàlues in (|2ü|l allows 
one to write the quantization condition (JTJ in the final form 

N 

r(n + 1 + í £ ) ((iV + l)/2) L /^ +1 ) 7fcL = 0. (23) 

L=0 

We have used the last expression of the quantization condition to find the lowest 
eigenenergies of the anharmonic oscillator (J2J for different vàlues of the coupling 
parameter g and four different choices of N. In the computation we have used a 
FORTRAN program with double precision. The results are shown in Tables 1 to 4. 

The procedure presented above assumes the capability to compute the N + 1 
coefhcients 7^ L (L — 0, 1, . . . , N) by summation of the series in (fTü|) . We have not yet 
proved rigorously the convergence of such series, albeit extensive numerical explorations 
guarantee its convergence for sufficiently large k, i. e., for n, in Eq. (|22j) . above a certain 
threshold which depends on the vàlues of the coupling parameter g and on the energy. 
Moreover, those explorations show that, the larger n is taken, the faster becomes the 
convergence. Investigations tending to elucidate that question are currently in progress. 

Besides the eigenenergies, our method determines also, in principle, the 
eigenfunctions. In the example considered, they are given by Eqs. (jHJ) and (|12|) . 
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Table 2. The four lowest eigenenergies of the Hamiltonian J2J|, for N = 5 and several 
vàlues of g. 



9 


Eo 




E-2 


Ep, 


-20 


- 11.56630147 


- 11.45854677 


0.56494700 


4.90729085 


-10 


- 2.83782675 


- 1.83075483 


4.90946147 


11.94279256 


-1 


1.03205834 


4.51533389 


10.48697985 


18.45464482 


-0.1 


1.27308185 


5.04058836 


11.08762465 


19.11537634 





1.29884370 


5.09787653 


11.15431820 


19.18880956 


0.1 


1.32441224 


5.15495387 


11.22099452 


19.26224408 


1 


1.54626351 


5.65933772 


11.81996788 


19.92310357 


10 


3.21711708 


9.93229322 


17.51589563 


26.43450876 


20 


4.48623513 


13.55329264 


22.99231828 


33.19354764 



Table 3. The four lowest eigenenergies of the Hamiltonian for N — 6 and several 
vàlues of g. 
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Ei 


E 2 


E3 


-20 


- 9.36607177 


- 9.13010587 


2.01035459 


7.97554684 


-10 


- 2.24187409 


- 0.87004433 


6.12159677 


14.16512836 


-1 


1.11369983 


4.84470202 


11.28130698 


19.99987959 


-0.1 


1.33949907 


5.33347217 


11.83181276 


20.59539382 





1.36377971 


5.38694202 


11.89300908 


20.66163760 


0.1 


1.38786579 


5.44024556 


11.95420520 


20.72789495 


1 


1.59799050 


5.91264617 


12.50470842 


21.32474109 


10 


3.22441873 


10.00630419 


17.83164730 


27.27876498 


20 


4.48680192 


13.57082013 


23.11371663 


33.63281210 



Table 4. The four lowest eigenenergies of the oscillator |J5J|, for N = 7 and several 
vàlues of g. 



9 


E 


Et 


E-2 


E;>, 


-20 


- 7.97489149 


- 7.59026706 


3.05916112 


10.19269195 


-10 


- 1.8474624 


- 0.17159144 


7.07320094 


15.87259291 


-1 


1.18393765 


5.12329191 


11.93911991 


21.26204013 


-0.1 


1.39832030 


5.58552094 


12.45475050 


21.81341553 





1.42143888 


5.63618503 


12.51210199 


21.87477520 


0.1 


1.44442247 


5.68671175 


12.56946066 


21.93615283 


1 


1.64542730 


6.13534277 


13.08581400 


22.48930458 


10 


3.23335919 


10.08415888 


18.13465608 


28.04433038 


20 


4.48835326 


13.59428939 


23.24781210 


34.07417453 
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Nevertheless, although the series in (|12|) converges for all finite x, it cannot be used 
safely for large vàlues of x, unless a considerable number of digits are maintained in the 
successive arithmetical operations. Certainly, the asymptotic expansion © can be used 
for sufficiently large vàlues of x (above about 5 units). But it is not clear the advantage of 
this procedure over the conventional numerical integration of the Schròdinger equation, 
especially if one needs the normalized wave function for a large number of points. 

To facilitate understanding the method, we have chosen above a very simple 
example: a one-dimensional anharmonic oscillator with only two terms in the potential. 
The procedure is equally applicable to isotropic D-dimensional oscillators with any 
number of integer powers of the radial variable in the potential and for any value of the 
D-dimensional "angular momentum" . It is also applicable, of course, to easier problems. 
Let us consider, for instance, three anharmonic oscillators algebraically solvable, namely, 
the Pòschl-Teller, the modified Pòschl-Teller and the Morse potentials. Their exact 
solution can be found in Ref. whose notation we adopt. In what follows, we 

concentrate on obtaining, by our procedure, the eigenenergies of the bound states. But, 
since the reflection and transmission coefficients are trivially related to the connection 
factors, our method is also useful for calculating phase shifts. 

In the case of the Pòschl-Teller potential 

, « 1 /k(«-1) à(à-1)\ Tr h 2 a 2 x , , 

V(x) = -Vo[ \, + \, , V = , k,A>1, (24) 



2 ysin (ax) cos 2 (aa:) J m 
defined in the interval x e [0, 7r/2], the Schròdinger equation can be written in the form 

. „ /l \ . 1 (k 2 k(k-1) A(A-1)\ , , 

in terms of a new variable 

y = sin 2 (o!a;) (26) 
and using, instead of the energy E, the parameter 

e = *£. (27) 

Equation (J25|) can be written as a hypergeometric one by means of the change of function 
done in Ref. Then, it is immediate to write the connection factors and to obtain 

the quantization condition. Nevertheless, let us ignore this fact and try to apply our 
method directly to Eq. í[25|). to be solved between the two regular singular points y = 
and y — 1. The solution regular at y — can be written as a power series 

oo 

u reg (y) = J2^y n+K/2 , a ^0, (28) 

n=0 

with coefficients given by the recurrence 



// [n—-+K 





(29) 
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The solution regular at y = 1 can be immediately written if one realizes that the 
differential equation (J23|) is invariant under the interchange 

í V \ í 1 -2/ 1 



and, therefore, 



k I A 



(1) Ü/) = EMl-í/) n+À/2 , &o/0, (30) 



n=0 



with coefficients given now by 



J n—1 




4a 2 



6n-2. (31) 



Both series in Eqs. (|28j) and (|30|) are convergent for y e (0, 1). We can, therefore, write 
a convergent (not merely formal, as in the problem considered before) expansion of the 
Wronskian. At the point y = 1/2, for instance, one has 

wr (Du Mo\ 1 ( f^an\ (m + X/2)b m \ 
W[u ieg , vS >] (y = 1/2) = -^-^ (í £ ^ J ^ E q ^=1 J 



+ (£ <^^(££,,. (32) 



y(o;) = -^ 7 a 2 __^ 2 . A>1. (34) 



\n=0 / \m=0 

Giving numerical vàlues to k and A, one can check that the Wronskian vanishes whenever 
k 2 

— = (« + A + 2n) 2 , n = 0,1,2,..., (33) 
cr 

as it should be. 

The modified Pòschl-Teller potential, defined for x G (— oo, +oo), reads 

h 2 a2 A(A - 1) 
2m cosh 2 (ax 
Instead of the (negative) energy E, we use the parameter 

^ - (35, 

Once again, the Schròdinger equation can be written as a hypergeometric one with 
adequate changes of variable and function. The change of variable used in Ref 
maps the interval (— oo, +oo) for the variable x onto [1, +oo) for the new variable. Here 
we prefer, however, to make a different change of variable, namely 

^ cosh 2 (ctx)' ^ 
in order to get the mentioned interval, where the differential equation has to be solved, 
mapped onto [0, 1]. The Schròdinger equation turns then into 



y\l - y)u" + y(l- |y) v! + I (-^ + A(A - l)y) u = 0. 



(37) 



Anharmonic oscillators 9 
The regular solution at the regular singular point y = can be written as a series 



; (y) = Y,*ny n+K/{2a) , a ^0, (38) 

n=0 

with coefncients obtainable by means of 

n („+í) c. = ((»-i+£) (»-5+s) - ^r 11 ) a »- < 39 > 

Now we need to write the well behaved solution at y — 1. Two independent expansions 
in power series of 1 — y of the form, 

oo 

u(y) = J2 h n (1 - y) n+ », Òot^O, (40) 

n=0 

with 

fi = or /i = 1/2 
and coeficients obeying the recurrence 

x> u fn, 1 , n2 , ^ 2 A(A-l) 



(ra+/x) (n+Ai-g) 6 «. = Í2 (n-l+^x) 2 + ^ - 4 - j 6n-i 



(n-2+/i) (n-|+Aí) - ^-il ] 6 n _ 2 , (41) 



2 

are physically acceptable. Those solutions with \i = and /U = 1/2 correspond, 
respectively, to even and odd wavefunctions in the variable x. Choosing the point 
y — 1/2 to evaluate the Wronskian of u reg and each one of those functions, one has 

ww t ,u](v = 1/2) = -^f{ (££)(£ Í!I ^=) 



n=0 / \m=0 

(42) 

\n=0 ^ / \m=0 ^ , 

It can be checked numerically that the right hand side of (J42j) becomes zero if, 

r „ , « „ K í A — 1 — 2n for even states, 

for n = 0, 1,2, . . . , < - = < A 

a [ X — 2 — 2n for odd states. 

We arrive finally to the case of the Morse potential 

V(r) = -D(exp(— 2ax) — 2 exp(— ax)), x=(r — r )/r , a > 0, (43) 



exactly solvable for angular momentum l = 0. By introducing, as in Ref. [Ï4] . a new 
variable 

27 

y = — exp(— ax), (44) 
a 



and denoting 



~ 2mErn 9 2mDri „ , , 

(3 2 = zr 1 , 7 2 = ^^, (3,1 >0, (45) 
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the Schròdinger equation becomes 

y 2 u" + yu' + [-4 + 2 2/ - 7^ J « = 0. (46) 

This equation presents a regular singular point at the origin and an irregular one at 
infinity. The physical solution, however, needs to be defined only between y = 0, 
corresponding to x — > 00 (r — > 00), and y = y = (27/a) exp(a), corresponding to 
x = — 1 (r = 0). Such physical solution must be regular at y = and become zero at 

y = y . The solution regular at y = can be given as a series 

00 

Uieg{y) = a nV n+l3/a i «0^0, (47) 

n=0 

with coefficients obeying 

n(ra + 2/3/a) a„ = -(7/a) a n _i + (1/4) a n _ 2 . (48) 

The other extreme of the interval of definition of the wave function, y = yo, is an 
ordinary point. There are two independent solutions of the differential equation, both 
finite at y — y^. But only the linear combination of them becoming zero at that point 
is physically acceptable. Let us call it Now, following our procedure, we should 
impose the cancellation of the Wronskian of M reg and w 1 ' at any point of [0, yo\. If we 
choose y = y , it becomes 

>VK cg ,« (1) ](!/ = 2/o)=M reg (!/o)« (1) '(2/o), (49) 
and, since u^' cannot vanish at y = y , the quantization condition reads 

Mrcg(Z/o) = 0, (50) 

an expression that could have been obtained trivially, without having recourse to our 
method. It is not difficult to see that, if one takes üq = 1 in (|47p. one has 

u Teg (y) = / /Q exp(-y/2) & ( í + — - 1 1 + 2^; y) , (51) 

V l a ca a 



and the quantization condition ()50|) coincides with that given in Ref. 

Unlike what happened in the case of potential (J2J), the numerical convergence of the 
power series giving the solutions of the Schròdinger equation in the three last examples 
is ràpid enough as to guarantee an accurate computation of the eigenfunctions. For 
instance, in the case of the Poschl-Teller potential, the series in (|28j) can be used for 
y < 1/2 and that in (JÏÏÜj) for y > 1/2, the coefficients a and b being determined by 
continuity at y — 1/2 and normalization in the interval y G [0, 1]. 
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